ASYMPTOTIC EXPANSIONS FOR ANISOTROPIC HEAT KERNELS 



LIVIU I. IGNAT AND ENRIQUE ZUAZUA 



Abstract. We obtain the asymptotic expansion of the solutions of some anisotropic heat 
equations when the initial data belong to polynomially weighted L p -spaces. We mainly 
address two model examples. In the first one, the diffusivity is of order two in some 
variables but higher in the other ones. In the second one we consider the heat equation 
on the Heisenberg group. 



1. Introduction 

This paper is devoted to the asymptotic behavior as time tends to infinity for the Cauchy 
problem associated with some anisotropic heat equations, with initial data in polynomially 
weighted L p spaces. 

There are a number of results on asymptotic behavior of the classical isotropic heat 
equation on the whole space: 

, ,t>o, 

(1.1) 

The key to obtain a complete asymptotic expansion of solutions as t — > oo is a decom- 
position of the initial data on the basis constituted by Dirac's delta and its derivatives, 
proved in [3] and that we recall for the sake of completeness. The asymptotic expansion 
reads as follows: 

Let G(-,t) be the heat kernel. For any 1 < p < N/(N — 1) and k > an integer the 
solution of problem (11.11) satisfies: 




\Lp(R n ) 



M*rdz)D°G(;t) <t-^||N fc+1 /ll 

rrf. al k Jr n j Lp ( rn ) 

\a\<k 

for any initial data f G L 1 (R N , 1 + \z\ k ) such that \z\ k+1 f G L P (R N ). 
We refer to [3] for a proof and various extensions. 

According to this expansion the asymptotic behavior of solutions of the heat equation 
can be described in terms of the Gaussian heat kernel and its space derivatives, using the 
moments of the initial data as coefficients. 

This decomposition is of isotropic nature; all space variables are treated equally, with the 
same homogeneity. This fact does not permit dealing with anisotropic heat kernels that 
would require different decomposition results, weighting the various euclidean variables 
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differently, or even with isotropic heat kernels, as the one associated with the heat equation 
(II. ip . but with initial data in anisotropically weighted spaces. 

In this paper we present new decomposition results for the initial data that lead both to 
new results for the heat equation ( II. ip and that, simultaneously, allow handling a number 
of anisotropic heat kernels. 

We focus mainly on two model examples to illustrate the key ideas and the phenomena 
that emerge due to anisotropy. 

The first one is the following diffusion equation, of mixed order, namely order two in 
some of the space variables, and order four in the other ones: 



;i-2) 



u t {z, t) = A 2 x u(z, t) + A y u(z, t), z = (x, y) G R m x R n , t > 0, 
u(z,0) = f(z), zeR m+n . 



The solution u can be written in convolution form as u(t) = Gt * f where Gt is the 
fundamental solution 



where the similarity profile G\ is a smooth function. 

Based on the decomposition obtained in [3] we get a new one (see Lemma 13.11 below) 
allowing us to prove the following asymptotic expansion. 

Theorem 1.1. Let k be a positive odd integer. For any f G L 1 (W n+n , 1 + |x| fc+1 + \y\i) 
the solution u of system ( 11. 2ft satisfies 



( 1 ^ I {p i) I r 

(1.3) u(-,t)- m .' ( / f&V^dxdy^D^G, 

<t-^\\(l + \x\ k+1 + \y\^)f\\ LHRN) 

Observe that the assumptions on the initial data weigh differently the various euclidean 
variables, taking into account the inhomogeneity of the diffusion operator. The resulting 
decomposition is in agreement with the different weight that the euclidean variables have 
in the fundamental solution, so that the overall contribution of its different derivatives is 
balanced. For simplicity we only have stated the result in the L 1 (M m+n )-norms but similar 
results can be stated in L p (W n+n ) spaces under some restrictions on the exponent p. The 
decomposition in [3] would be insufficient to obtain such results since it is of isotropic 
nature, thus weighting equally all space variables, which is not sufficient to deal with the 
anisotropy of the heat kernel under consideration. The method presented in this paper 
allows treating separately the euclidean variables with different homogeneity. 

As we shall see, the decomposition used to obtain this result can also be applied in 
the case of the classical isotropic heat equation, leading to new results on its asymptotic 
decomposition as t — > oo. 
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The second problem we address is the heat equation on the Heisenberg group: 

u t {z, 9, t) = A m u(z, 9, t) (z, 9) G R 2n xl,{>0, 

u(z,9,0) = f(z,6), , (z, 9) G M. 2n x R. 

Let us recall that, as for the classical heat equation, the heat semi-group on HP is given 
by the convolution between the initial datum and the fundamental solution, a function in 
the Schwartz class (see jl]). 

Theorem 1.2. (T4jj There exists a function H G 5(lR 2n+1 ) such that, u, the solution of 
the heat equation on the Heisenberg group (1 1.4ft is given by 

U (;t) = f*H t , 

where * is the convolution on the Heisenberg group while H t is defined by 

Moreover, the function H can be computed explicitly (cf. [Ijj 

H{{z,9)) = -^— I ( w T exp(^-I^-^- W 
VV ' ;; (47r)™ +1 J R Vsinh(2cr)/ P V 2 2 tanh2a/ 

In this case we obtain the following asymptotic expansion. 

Theorem 1.3. For any f G L 1 (IR 2n+1 , 1 + \z\ 2 + \9\) the solution u of problem ( 11.41) with 
initial data f satisfies the following 

(1.5) u(-,t)-[ I f)H t + y j [ I z j f(z,9))Z j H t 



<rl(i + |z| 2 + |0|)/|| L i (R2n+1) , t>o, 

where Zj, j — 1, . . . , In are the first 2n vector fields of the Lie algebra on HP. 

The results presented here can be extended to more general situations, for example to the 
heat equation on homogeneous Carnot groups. We refer to [6] where this issue is addressed, 
although the results in that paper do not take into account the different homogeneity in 
what concerns the needed assumptions on the initial data. Note also the existing results 
for heat equations in heterogeneous media. For instance in [5] by combining the results 
in [3] and Bloch-wave expansions a complete asymptotic expansion is given for the heat 
equation in periodic media. We also refer to [2] for similar expansions in the context of 
Vlasov - Poisson - Fokker - Planck equations. 

The paper is organized as follows. In Section [2] we consider the classical isotropic heat 
equation to illustrate that the use of the new decomposition results we present here leads 
to new asymptotic expansions as t — )• oo that can not be derived out of [3] . In Section [3] 
we consider the anisotropic model (jl.2p . proving Theorem 1 1 . 1 1 and we compare it with the 
asymptotic decomposition results one could get out of the isotropic decomposition results 
in [3]. In the last section we present some classical facts about the Heisenberg group and 
give the proof of the corresponding asymptotic expansion stated in Theorem 11.31 
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2. Decomposition Lemmas 

In this section we obtain some variants of the results in [3] on the decomposition of 
functions defined in the euclidean space in the basis of Dirac's deltas and its derivatives. 
We shall derive some applications to the asymptotic behavior of the isotropic heat equation 
for initial data in anisotropically weighted spaces. 

In what follows we denote by X>(R ) the space of C°° and compactly supported functions 
and by V'(R N ) its dual. 

We first recall the following decomposition [3]. 

Lemma 2.1. (\3\) Let k be a positive integer. For any f E V{R N ) it holds 

(2.6) f(x)x a dx)D a 6 + V D a F a , 

where the functions F a , defined for |a| = k + \, are given by 

(2.7) F „ (l) = (i + 1) (z^l J f 1 (1 _«).(£)./(£)* 
Moreover, for any p E [1, jfzj) the following holds 



^r + 7iUp(R^)- 

Remark 1. By density, the above result holds for all f E L 1 (M. N , 1 + \x\ k ) with |x| fc+1 / E 
L p (R N ). 

The above result is a consequence of Taylor's formula with integral remainder, applied 
to any test function ip: 



(2.8) <p{z) = T (D*<p)(0)£ + (k + l) ~y A 1 - t) h (D"<p)(tz)dt, z E R N . 

al al / n 

|a|<fe |a|=fc+l u 

We now derive some anisotropic variants. For, given z E ~R N , N > 2 we decompose it as 
z = (x, y) E W 71 x IR n with m + n = N. We then use a splitting argument and apply first 
Taylor's expansion in the y variable with x fixed and then perform Taylor's expansion in 
x. 

Lemma 2.2. For any (p E T>(R N ) the following holds 

^) = E(£»(0)^+E [(fc + l-|7l) E J\ j\l-t) k+bl (D^)(tx,0)dt 



\a\<k \<y\<k |/3|=fc+l- 7 

,7 r 1 



+ (fc + l) V 4 / z=(x,y)ER N . 



|7|=*>+1 
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Proof. Using Taylor's expansion in y with x fixed we get 

<p(z) =<p(x,y) = 0)^ + (* + l) ^ f\l-t) k (DJ<p)(x,ty)dt. 

h\<k 7 ' \y\=k+l 7 ' Jo 

Now, applying Taylor's expansion with respect to x we get: 



/3! 7! 

| 7 |<fc |/3|<fc-|7| 



| 7 |<fc |/3|=fc+l- 7 P ' J ° 

+ + E ^j\l-t) k (D^)(x,ty)dt. 



h\=k+l 

Writing a = (^,7) we obtain the desired identity. □ 

We now introduce the space L p (R m , \x\W] L 1 (R", |y| l71 )), constituted by the functions / 
such that \x\W\\\y\^f(x r )\\ L i {Rn) belongs to L p (R m ), i.e.: 

( 2 -9) ll/H^PCMm .l r ||/3|. ririRn I„|l7h1 = / M^'ll/O^) OllrlnRn-U/llTh^ < °°- 



LP(R m ,|x|l^l; L 1 (IR™,|»/|ItI)) / 1^1 H ^ ^' ' H L 1 (R";|?,|H)'- 

Following the same ideas of treating differently the different variables x and y we get 
the following result. 

Lemma 2.3. For any / e V(R m+n ) it holds 

(2.10) f(x,y)= ( "i;' ( ;' 7) ' ( / /(*, ^V^dy)^ 



+ E E ^^^/W^Wi/) 



| 7 |<fc |/3|=fc+l-| 7 | 7 ' 

+ ]T D^F 1 (x, y), (x,y)eR m xW\ 

h\=k+l 

where the remainder terms are defined as follows: 

^ ) = (z3» J f 1 (1 -«)«-'(f)'X./(f,^ 

/or ItI <A;,|/3| + |7| =k + l, and 

F 7 (x, y) = (k + 1)^^ ^(1 _ t) fc (|)V(x, f )* 

/or |7| = k + 1. 
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Moreover, for 1 < p < it holds 

(2-H) ||F 7 || LP(]R ]V) < +1 /||iP(RJV) 

and for 1 < p < it holds 

» — ^ m— 1 

(2-12) ll[^/WU p (R m ) - ll/ILp(R™,|a;|^l; L 1 (IR n ,|j/|lTl)) ■ 



Remark 2. Tae mam term is the same as in Lemma \2. 1\ since 

(-l)'^ 7)l f f f^y )x ^ dxdy ) D ^ So= J2 [ f(z)z a dz)D a 5 . 

\W)\<* P1 JRN WRiv 

estimating the remainder terms _F 7 and [M/]^ different anisotropic spaces from those 
of Lemma \2.1\ are needed. 

By density, the decomposition holds for all f G L 1 (M. N , 1 + | • | fc ) that belongs to 

XP= p| L p {R m , |x| fc+1 ~' 7 '; L 1 {R n ,\y\^))[]L p {W n - L p (M n , \y\ k+1 )). 

h\<k 

We also define the natural norm on X v 

ll/IU p := E ll/IL§(K' n ,|a;| fc + 1 -lTl; £,l(K»,|y|H)) + || M^Vlll^RN). 
0<|7|<fc 

Whenp = 1 the space X p is exactly the one needed in Lemma \2. 1\ i.e. L l (M. m+n , \ (x,y)\ k+1 ) . 

Proof of Lemma\EM We write R N as l m xl n . We first fix x G R rn and apply Lemma O 
to the function y — > f(x,y). Then 

(2.13) /(^2/)=r^-(/ /(*, v)v 7 dvW*o(v) + E D1 v F 7( x >y)> 

h\<k T Wr " | 7 |=*+1 

where 

(2-14) F 7 (x,n) = (fc + l)til^ [\l-t) k ( y -yf(x, y -)^. 

1- Jo 1 11 

A similar argument as in [3] shows that for any x G M m and p G [1, -^j) we have 

(2-15) 11^(^0 ||£P(R») < III ■ | fe+1 /(^,-)IU,(«"). 

Taking the L p (R m )-norm in the x- variable we obtain estimate ( 12. lip . 
We denote by [M/] 7 (x) the y- moment of order 7 of f{x,-) in R n : 

(2.16) [Mf] 1 (x)= [ f(x,y)y^dy. 



We apply Lemma 12. II to the function x — > [Mf] 1 {x) and obtain 

[MfUx)= Yl [M/] 7 (x)x 8 dx)L»f5 (x)+ £ 

\p\<k— \j\ |/3|=fc+l-|7 



/3 7 



;r 
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where [R-f]p y is given by 

Writing explicitly the function [Mf] 7 , we have 
(2.17) 

[MfUx)= f{^yVy 1 dxdy)D%{x)+ Yl ^PVW*) 

\l3\<k-\l\ P ' JRN \p\=k+l-fr\ 

where [R-f]^ are given by 

Replacing [M/] 7 in ( 12 .131) with identity (12.17ft we obtain 

f = H E \ m [ / f(x,y)x^dxdy)D%(x)D;S (y) 

+ E E ^r-^[^/]^(^)^o(y) + E ^W^v) 

| 7 |<fe |/3|=fc+l-| 7 | ™" | 7 |=fc+l 
-1)1(^,7) 

l(0.7)l<* 

+ E E ^^DiiKfifhWW + E W^)- 

| 7 |<fc |j8|=fc+l-[7| M=fc+1 

Proceeding as in [3] and using Minkowski's inequality, for any 1 < p < we get estimate 
(12.121) on the remainder [7£/]/3 7 



||[^/]^ 7 ||lp(r"») < MM f(%,y)y J dy\ 

The proof is now finished. □ 



We now apply the new decomposition obtained in Lemma I2TB1 to the heat equation (11.11) . 

Theorem 2.1. Assume that N = m + n. For any 1 < p < min{^-j-, ^ry} and k > 
integer, the solution of problem ( 11. ip satisfies: 



"M - E / /W^)^(-,t) 



|a|<fe 

/or any zratoaZ data / G L 1 (R JV , 1 + |z| fe ) n X' p . 



Lp(R n ) 
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Remark 3. In the case p = 1 the above estimate is the same as the one obtained in [3] by 
applying Lemma \2.1\ 

The range of exponents p for which the above Theorem applies is larger than that in [3] 
but the space X v where the initial data is taken is more complex and of anisotropic nature. 

Proof. We now apply Lemma [2.31 to the initial data /. Then 



G(t)*f - J2 



\a\<k 



f(x)x a dx)D a G(t) 



LP 



^ E E || G '( t )*(^ [7e/] ^ (a;)jD ^ o(2/) )IL,(R-) + E \\G(t)*D2F 7 (x,v)\\„ 



0<|7|<fc \P\=k+l-\y\ 

s E 

0<h\<k 

* E 

0<h\<k |/3|=fe+l-]7| 

s E 

0<|7|<fc 

s E 

0<|7|<fc 
fc+i 

< t 2 || / \\xv. 



|7l=fe+l 



\\DiDZG{t,',v)* a [Kf\ fh {-)\\ l , 0l » ) + \\DJG(t)*F 7 (x,y)\\ LP 

=fe+l-|7| | 7 |=fc+l 



E |ii^ 7G wiUiii^/]^iiLg 



+ t -— R 



7 \\LP 



E ll^/WWI^^t)^ + \t\*\\\y\ k+1 f\\L* 

=fe+l-|7l 

E i< 

=fc+l-| 7 | 



2 ^ pJ ||/||L P (ffim H |^|. Ll(R»,| y |l7l)) + \t\ 2 ||M + /||lp 



The proof of Theorem 12.11 is now finished. 

3. An anisotropic heat equation of mixed order 
This section is devoted to the analysis of the following system 

u t (z, t) = A 2 x u(z, t) + A y u(z, t), z = (x, y) G R m x W 1 , t > 0, 



□ 



(3.18 



u(z,0) = f(z), 



z e 



ira+ti 



The solution u of the above system is given by u(t) = G t * f where G t is the fundamental 
solution given by 

Gt (x,y)=r?-5G 1 (^ i ,^) ) 
and the self-similar profile G\ is so that its Fourier transform is such that 

G 1 (e,r 7 ) = e-I«l 4 -H 2 . 

We now apply Lemma [27T1 to the initial data /. For any 1 < p < N/ (N — 1), N = m + n, 
we obtain that u, solution of (I3.18p . satisfies 

(-l)H 



LP 



f\\LP(R N ) E ^"^tlU 1 



|a|=fc+l 
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or, writing each index a of length N = n + m&sa = (/?, 7) G Z m x Z n , 



*)- E 



I09,7)l<* 



(_l)l(0.7)l 

^7! 



f(x,y)x^dxdy)D^G 



LP 



<\\z\ k+1 f\\LP(K») E n^ G *iU ] 

|/3| + | 7 |=fc+l 

Remark that, due to the anisotropy of the fundamental solution, the derivatives in the 
x- variable decay differently from those in the y- variable: 



(3.19) \\D^G t \\ Lq{ 
As a consequence, for large t we have 



J£l_M_iV/i_i 



£ 4 2 2 



a-i) 



? G [l,oo]. 



■*)" E , ( / J{x,v)xPyidxdy)D%G t 

1 09,7)1 <* 

<ikl" +1 /llL, 



LP 



E 

|/3| + | 7 |=fc+l 



,_L£i_N ^ . fc+i ,| 1 ,l. , 1 



LP 



Note however that in the left hand side of the above inequality there are terms that decay 



faster than the remainder t 4 . Accordingly it is natural to split the terms in the left 

fc + 1 

hand side in two parts: those that decay faster and, respectively, slower than t ~. 
Define 

A(p, k) = {(a, b) G I? : a, b > 0, k + 1 - 2JV(1 - 1/p) < a + 2b, a + b < k}. 

It follows that, as long as / belongs to some weighted spaces, the solution u of system 
(13. 18)) satisfies for large time t the following estimate 



(3.20) «(-,*)- E 

|/3|+2| 7 |<fc+l-27V(l-l/p) 

<r^\\\z\^f\\ LP 



(_1)I09.7)I 
^7! 



f{x,y)x^dxdy)D^G 



LP 



<t* 



||; 2 |fe 



£ 4 "a" Tv p- 1 



I LP 



+ E 

C8,7)eA(p,fc) 

+ E 117 lLJ,l - : % -l-H " l.v ' 1 
08, 7 )eA(p,fc) 



L 1 (R JV ,k|l /3 l|^l l71 ) 



We now state rigorously the result to which this analysis leads. For the sake of simplicity 
we state it for p = 1 only. 
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Theorem 3.1. For any f G L 1 (M. n+m , 1 + \x\ k+1 + \y\ k+1 ) the following holds for large time 



u{ 



■«)- E 



|/3|+2| 7 |<fc 



-1)108,7)1 



f(x,y)x^dxdy)D^G 



Proof. We apply estimate (13.201) to obtain the following: 

'_l)l(/3,7) 



/3! 7 ! 



|/3|+2| 7 |<fc 
_fc+i 



f(x,y)x^dxdy)D^G 



+ E 

(/3, 7 )eA(p,fe) 



,|x|l/%|M) 



fc+i 



< r^- + lyO/lkw + ||(N fe + M fe )/|| L1 

Since / belongs to L 1 (M Ar ) we obtain the desired result. 



□ 



Observe that, in view of (13.191) . we need twice more derivatives in x than in y to obtain 
the same decay of the kernel G t - We now give a new decomposition that, when applied to 
the anisotropic equation (I3.18p . distinguishes the two variables x and y. 

We consider the case when k is an odd integer, k + 1 — 2m + 2. Since we have seen that 
\\DxG t || lp (r n ) — \\D^G t \\ LP ( R N-), \a\ = 2|/3| = 2m, we will take into account that, for the 
decay rate of the kernel, one derivative in y has the same effect as two derivatives in x. 

Lemma 3.1. Let k be a positive odd integer. For any f e V(W m+n ) it holds 
(3.21) 

(_1)I09.7)I 



f(x,y)= 

|/3|+2| 7 |<fc 



/3! 7 ! 



+ 



E 



E 



f(x,y)x l3 y' y dxdy)D^5o 

:_i)H 



7! 



Dl [1lf] M {x)D^ {y) + D V F ^V) 



|7l<(fc~l)/2 |/3|+2| 7 |=fc+l 

where the remainder terms are defined as follows: 



| 7 |=(fc+l)/2 



fc + l(_lff h 
2 7! 



t 



t't r ' 



for |7| 



fc+i 



, and 



(_!)!/»! 



defined for | 7 | < (Jfe - l)/2, |/3| + 2| 7 | =Jfe + l. 
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Moreover, for 1 < p < it holds 



k+l 



(3.22) \\F,\\ LP{ m) < \\\y\— f(x,-)\\LP{^), 

and, for 1 < p < 

(3-23) ||[^/]^ 7 |Up(IR™) < ||/||lp(R™ ) [x[IPI; Ll(R",|2/|W))- 

Remark 4. i?y density the results hold true for a larger class of functions f . 

Proof of Lemma \3. 11 Set k + 1 = 21 + 2. Fix x £ M m and apply Lemma I2TT1 to the function 
V -»■ /(»,!/)■ Then 



(3.24) /( x , y ) = ^LiffC / /(x,y)^)^ (y)+ V) D^ y F 7 (x,y), 

h\<i T ^ l7l=m 



where 

T) m Z" 1 /, y,dt 



(3.25) F 7 (x, y) = (J + 1)^— J (1 - (f) 7 /(*, |)^, | 7 | = * + 1- 
Observe that for any x £ IR m and p £ [1, -^-) we have 

(3.26) \\F,(x,-)\\ LP(m < \\\y\ l+1 f(x,-)\\ LPm . 

Taking the L p (R m )-norm in the x variable we obtain estimate (13.221) . 

Recall the definition of [M/] 7 given in (12.161) . We apply Lemma 12.11 to the function 
x — > [Mf]y(x), with a remainder of order fc + 1 — | 7 | and obtain that 
(3.27) 

WfUx)= i - Z W~{[ [Mf\ 1 {x)x?dx)D%{x)+ °i[Kf\fh( s 

|/3|<fc-2| 7 | ' jRm ' |/9|=fc+l-2|7| 

where 



[X 
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Using decomposition (I3.27j) in (I3.24p we obtain that / can be written as 



(/ f(x,y)x^dxdy)D%(x)+ 



h\<m ' |/3|<fe-2|7| K JK 

+ Yl DP[Kf]p 7 (x))D15 (y) 

|/3|=fc+l-2| 7 | 

+ J2 D v F ^y) 

| 7 |=m+l 

_1)M 



E E 

|7|<m |^|<fc-2|7| 



7! /3! 



/(*, y)x^dxdy D%(x)DMy) 



+ E E 

|7|<m |/3|=fc+l-2| 7 

(_1)I(A7)I 



I) 1 



0"! 



-Dg[Kf\ fh (x)DZ6 {y)+ ^ ^(ar.y) 



E 



|/3|+2| 7 |<2m+l 



/3! 7 ! 



| 7 |=m+l 

f(x,y)x^dxdy)D^5 



+ E E ^^ft^/w^w + E D v F ^y)- 



\-y\<m |/3|=2m+2-| 7 | 



7 



| 7 |=m+l 



Estimate (13.231) on the remainder [JZf]^ is obtained as follows: for any 1 < p < we 
have 



ll^/KIU^R™) < W 9 / f(x,y)y J dy 



< 



Lp(R") 



LP(R™,|z|l' 3 l; L 1 (R m ,|3/|M))- 



The proof is now complete. 



□ 



We now apply this new decomposition to obtain an asymptotic expansion for the solu- 
tions of (13.181) . and prove Theorem 11.11 
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Proof of Theorem ! We now apply Lemma [3TT1 to our initial data / G L 1 (lR Ar , l + |a;| fc+1 + 
Then 

G9.7) 



E J, ( / f(x)x f3 y^dxdy)D^ y G 

E E (^/W^W 



I/ 1 (R JV ) 
L 1 (IR JV ) 



|/?|+2| 7 |<fc 

< 

| 7 |<(fc-l)/2 |,8|+2|7|=fc+l 

+ ]T ||G**^F 7 (a;,y)|| LW 

| 7 |=(fc+l)/2 

< E E \\D^ y G t (-,y)* x [nf]^(-)\\ LHRN) 

l7l<(fe-l)/2 |^|+2|7|=fe+l 

+ \\D;G t *F,(x,y)\\ LHm) 

| 7 |=(fe+l)/2 

< E E ||p^ii Li(Rm) n[^], 7 iu i(Mm) E rai 



r 1 /]fgn A 

0<|7|<& |^|=fc+l-|7| ' | 7 |=(fc+l)/2 

< E E ii m^hi\\D^ y G t \\ LHRN) + itr 4 * 1 iiM 4 * 1 /iu W) 

|7|<(fc-l)/2 |/S|+2| 7 |=fc+l 

V-^ 1 I l_J£l±£Mll nil I I fc+J- ni i fc+1 -I, 

< 2^ E 1*1 4 II/IUiokMx|I/»imN) + 1*1 4 II Is/I 2 JIUhk^) 

l7l<(fc-l)/2 |/3|+2| 7 |=fc+l 

^r^lKi + N^ + lyl 4 ^)/!!^,. 

The proof of Theorem 11.11 is now complete. □ 

Remark 5. Remark that all the terms in the left hand side of the main estimate of Theorem 
11.11 have a slower decay than those of the right hand side. Also observe that, with respect 
to Theorem \3.1[ in Theorem \l.l\ we have assumed only integrability of (k + l)/2 moments 
in the variable y. 

4. The Heisenberg Group 

First of all we give the definition and describe the main aspects of the Heisenberg group. 
Given (z, 9) of M 2n+1 in the form (z, 9) = (z±, . . . , z n , z n+1 , . . . , z 2n , 9), the composition law 
on the group is given by 

n 

(z, 9) o (z 1 , 9') = ( z + z ',9 + 9' + 2 5>„ + i^ - 

3=1 

In this way (R 2n+1 ,o) is a Lie group, whose identity element is the origin (z, 9) = (0,0), 
the inverse being given by (z, 9)^ = (—z, —9). 
Let us now consider the dilations 

5 X : R 2n+1 -> R 2n+1 , S x (z,9) = (\z, \ 2 9). 
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We have that 5\ is an automorphism on (IR 2n+1 ,o) for every A > 0. In this way HP = 
(R 2n+1 , o, 5\) is a homogeneous space. The Lie algebra on HP is given by the left invariant 
vector fields 

Zj = d Zj + 2z n+j d e , Z n+j = d Zn+j - 2zjd e , j = 1, . . . , n, 6 = d e . 
The convolution product of two functions / and g on HP is defined by 

{f*9)(w)= / f(w o i r l )g{y)dv = / f(v)g(v~ 1 o W )dw, 



where v _1 is the inverse of element v with respect to the group operation o on HP. Here, 
dw is the Haar measure on HP which is exactly the euclidian measure on R 2n+1 . It should 
be emphasized that the convolution on the Heisenberg group is not commutative. However 
if P is a left invariant vector field on HP then one sees easily that 

P(h*h) = h*Ph = Ph*f2- 

The canonical sub-laplacian on HP is given by 

2n 

In the case of the Heisenberg group HP the exponential and logarithmic maps (see [1] , Ch. 
1, p. 48, for the definition of the Exp and Log maps on general homogeneous Lie groups) 
are given by (see pQ, p. 167 and Remark 3.2.4 on p. 163) 

In 

(4.28) Exp ( z i Z i + 9Q ) = (M), * = ■ ■ ■ , 

3=1 

and 

2n 

(4.29) Log(^, 0) = J2 z i Z i + 9Q - 

3=1 

We point out that properties (I4.28P and (I4.29P are not true in general for homogeneous 
Lie groups as we can see in the following example. Following [lj, Section 3.2, p. 158 let 
us set R N = R m x ]R n and denote its points z = (x, 6) with x G lR m and 6 e lR n . Given a 
m x m matrix B with real entries let 

(x,6)o(Z,t) = (x + Z,6 + t + ±{Bx,0), 

where (•, •) stands for the inner product in R m . Again 5\ : R N — > R N , 5 x (x, 6) = (Ax, A 2 ^) 
is an automorphism of (R N , o) for any A > 0. Then (R , o, 8\) is a homogeneous Lie group 
of step two. In this case the Exp and Log functions are given by (see (3.11), p. 167 in PQ) 

Exp((£,t)-Z) = (£,t + -(££,£)) 



and 



where 
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(£,r)-z = E^ + E T ^ 

j=l i=l 

and Xi, . . . , X m , Q%, . . . , 9 n is the Jacobian basis. 

We can see that (I4.28[) and (I4.29P hold if and only if B is a skew-adjoint matrix. In the 
case of H n this holds since the matrix B is given by 



B = 4 



In 
-In 



In the particular case of skew-symmetric matrices B we believe that the results of this 
section are still valid. 

We now recall some results concerning Taylor's formula with integral remainder on ho- 
mogenous Carnot groups (G, o). 

Lemma 4.1 (Lemma 20.3.8, pQ). Let x, g G G and u G C m (G, R). Then we have 

m „i 
(4.30) u(xoh) = V-((Log/i) fc M)(x) + — / (i-s) m ((Log/i) m+1 w)(xoExp(sLog/i))c2s. 

k=0 J0 

Let us now apply this Lemma to the case of the Heisenberg group H n . Using the structure 
on H ra we obtain that for any h G H n and sGl the following holds 

Exp(sLog h) = sh. 

Thus Lemma [4. II (see Corollary 20.3. 9 in [1J, p. 755) gives us the following: 
Lemma 4.2. Let x,/i6l" and u G C m+1 (H n , R). Then we have 

m 2n+l , 

K_ 

k\ 



u(xoh)=J2 E ^r^(Z tl . . . Z lk u)(x) 



k=0 ii,...,ife=l 

2™+! ;. U rl 



+ E ferxrTT / ( x - s ) m ( z « • • • ° sW*, 

ii £ ^ ^° 

tl,.. - ,t m -|-l — 1 

where h,2n+i = $ an d Z 2n +i = ©■ 

We now follow the same ideas as in the case of the anisotropic heat equation making 
partial decompositions in some of the involved variables. Here, in the context of the 
Heisenbergh group, variables z = (z\, . . . , £ 2n ) and z 2n +i = 9 will be treated separately. 
For simplicity we treat here only the case m = 1. 
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Lemma 4.3. For any if G I>(R 2n+1 ) the following identity holds: 

2n 

<p(z,0) =^(0,0) + ^z,(Z^)(0,0) 



3=1 

2n 



«1 m „x 

+ / (1-s) V) z j z k (Z j Z k <p)(sz,0)ds+ / 0(e<^)((z,s0))cis. 
Remark 6. We emphasize the difference with formula (20.63) in [I], p. 761: 

2n 2n „\ 

<p{z,0)=<p(O,O)+J2 z j( Z 3 ( P)( o ,°) + d ( e( P)( o > o )+ / 0--s)zjZ k {ZjZ k <p)(sz,se)d8. 

3=1 j,k=l Jo 

Proof. We apply Lemma [4.21 to x = (z, 0) and h = (0, 9), elements from R 2n x R. Thus 

(p(z,6) = (p(z,0)+ [ (0ep)((z,O)o(O,s9))ds = v(z,O)+ I 6(Q(p)((z,s6))ds. 
Jo Jo 

We now apply Lemma I4T21 with x = (0, 0) and h = (z, 0): 

2n „i 2n 

<p{z,0) = <p{0,Q)+ y %2z j {Z j <p){0,0)+ / {l-s) y £z j z k (Z j Z k <p)(sz,0)ds. 

3=1 J ° j,k=l 

Putting toghether the above identities we obtain that 

2n „\ 2n 

(p(z,t) =<p(0,0) +J2 z j( Z 3<f)(0,0) + / (1 -s) z j z k (Z ] Z kV )(sz,0)ds 

3=1 Jo j,k=l 



+ / e(e<p)((z,se))da 

Jo 

and the proof is now finished. □ 

We now apply the Taylor like identity in Lemma 14.31 to obtain a decomposition of a 
function / into Dirac's delta basis adapted to the Lie algebra on HP. 

Theorem 4.1. For any f e L x (R 2n , 1 + \z\ 2 + \9\) the following decomposition holds: 
(4.31) 



2n „ 2n 



f = ( I /Wi] ( / z j f(z,e)dzd8)z j 5 + (QF) + f] Z 3 Z k (F 3 M9)) 



w/iere 

(4.32) F(*,0) = - [ -f(z,-)ds 



s s 



and 

, i „ . s f f 1 , , Zj z k n ,z 9.dsd9 

(4.33) fy(*) = / / (l- s )-l^/(- -)_. 



o 
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Moreover, 

||-nUi(R2»+i) < ||^/||li(ir2u+i) 

and 

||-^jfc||L 1 (IR 2n + 1 ) ^ lki 2 fc/IU 1 (R 2 "+ 1 )- 

Proof. First of all it is sufficient to prove the above theorem for smooth / and then extend 
the result for any / G L 1 (R 2n+1 , \z\ 2 + \9\). Let us choose / and if two smooth functions. 
Applying Taylor's expansion obtained in Lemma [4.31 we get the following: 

» 2n 

(f^)v,v= / f(z,e)L(o,o) + J2 Zj (z j(p )(o,o))dzde 

« «i 2n 

+ 



+ 



p\ Zn 

f(z,9) / (1 — s) y ZjZ k (ZjZ k cp)(sz,0)dsdzd9 

[ f(z,6) [ 0{Qip){z,s6)dsdzd8 
Jm 2 ^ 1 Jo 

< (5(o,o), f > / 70, 0)<fed0 - y~] < ZjSfpfl), f> f(z, 6)zjdzd0 

+ f; / (z jZkV )(z,o) [\i-s)f( z -,ef-i z ^dzde 

~ k ~Z 1 JR 2n + 1 Jo s s s s 

+ [ [e<p)(z,0) [ -f(z,-)-dzd9. 

,/ R 2n + l Jq S S S 



Observe now that 



{ZjZwKz, 0) f\l - *)/(-, QfJ-^^dzdQ 
Jo s s s s 

jR 2n l Jm.Jo s 



s s s 2n 



dz 



(ZjZ k (p)(z,0)F jk (z)dz = (F jk 5 (9), ZjZ k tp) v > iV 
(ZjZ k (F jk 5 o (0)) , (p)w,T> 



where 



Denoting 



F jk {z) = [ [ ( 
Jr Jo 



l-s)f(-,9) 



z n ^Zj z k dsd9 
' s s s s 2n 



n:..H) j -f(z,-)ds 



-1 s s 

we obtain the desired decomposition. The estimates of the L 1 -norms of F and Fj k imme- 
diately follow. □ 
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We now apply the decomposition obtained in Theorem 14.11 to obtain the asymptotic 
behaviour of the solutions of the heat equation on the Heisenbergh group HP. 

Proof of Theorem \1.3[ Using decomposition (I4.3ip we have that 



u{z, 



e,t)-([ + zjffrejjZjHt 

\J R 2 n + l / " \J R 2n + l / 



2n 

< ||6F * iJi|| L l( R 2n+l) + \\Z j Z k {F jk 8o{9)) *i/ t || i l( K 2n+l) 

j,k=l 

2n 

= \\F * 6#t||l,l(R2n+l) + ||(F ifc 5 O (0)) * (ZyZfe^Hil^n+l) 

j,fc=l 

2n 

< ||-^||l 1 (K 2 «+ 1 )||@-^Jl 1 (IR 2 ' 1 + 1 ) + ^ ||^fc||Ll(K 2 ")ll'^i^^lU 1 ( R2 ™ +1 ) 

j,k=l 

2n 

< ||^/|Ul(K2n+l) + 2^ lkj 2; jfc/||Li(R 2n + 1 ) - 11(1^1 + k| 2 )/l|L 1 (K 2n + 1 )- 

The proof of the decay rate property is now finished. □ 
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